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Monotonicity  Conditions  of  Curvature 
for  Bezier-de  Casteljau  Curves 


Jean-Charles  Fiorot  and  Laurent  Schiavon 


Abstract.  In  this  paper,  we  deal  with  the  monotonicity  of  curvature 
problem  for  Bezier-de  Casteljau  curves.  We  focus  more  particularly  on 
the  cubic  case.  A condition  about  decreasing  curvature  at  the  origin  of  a 
cubic  curve  is  given  so  that  it  implies  decreasing  curvature  at  every  point. 
The  corresponding  cubics  are  determined  by  their  control  polygon. 


§1.  Introduction 

The  problem  discussed  here  concerns  the  shape  control  of  curves;  mainly  how 
to  obtain  curves  with  a monotone  variation  of  curvature.  The  aim  is  to  find 
the  widest  class  of  curves  with  monotonely  increasing  or  decreasing  curvature 
variation.  This  problem  mainly  arises  in  car  body  shape  design. 

Let  P be  a n-degree  Bezier-de  Casteljau  planar  curve  defined  on  [0, 1] 
given  by  its  control  polygon  {Po>  Pi,  ... , Pn}.  Let  us  consider  l =|  PqPi  |,  and, 
Vi  = 1, ..  .,n-l, set  ft*  =|  P,Pi+i  | / | Pi- ip  |and<p,  = (PjPj+i, Pj_iPj).  We 
say  that  P admits  the  representation  (hi, . . . , hn-i;  'Pi,  • • • , Pn- i?  /)•  Higashi, 
Kaneko  and  Hosaka  [3]  characterized  the  monotonicity  of  curvature  when  Vk  = 
1, . . . ,n  — 1,  hk  = h and  ipk  = <p  '•  h cos  p > 1 and  h < cos  ip  are  respectively 
the  condition  of  decrease  and  increase  of  the  curvature.  This  model  has  been 
used  recently  by  Mineur,  Lychah,  Castelain  and  Giaume  [4]  to  control  a shape 
when  fitting  a curve  to  a set  of  given  data. 

Here,  we  focus  on  the  cubic  case.  From  the  representation  (hi,  /12;  Pu  P2] 
l)  of  P,  we  apply  the  de  Casteljau  Algorithm  at  t value  belonging  to  the  interval 
of  definition,  and  we  determine  by  induction  the  parameters  characterizing  the 
two  segments  given  by  the  Subdivision  Algorithm.  This  process  enables  us  to 
determine  the  curvature  p.  Then,  we  seek  cubic  curves  for  which  decreasing 
curvature  at  the  origin  implies  decreasing  curvature  at  every  point.  Such 
curves  are  determined  via  the  parameters  r*  = hi! hi,  <pi  and  <p2-  The  study 
falls  into  two  cases  : r*  > 1 and  r*  < 1.  For  the  second  case,  we  give  a more 
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strict  decreasing  curvature  condition  at  the  origin  in  order  to  get  the  decrease 
of  p everywhere. 

In  the  framework  of  Fiorot  and  Jeannin  [1]  and  Fiorot,  Jeannin  and  Taleb 
[2],  an  attempt  to  extend  these  results  to  rational  cubic  curves  has  been  made. 


§2.  de  Casteljau  Algorithm 

Let  P be  a Bezier-de  Casteljau  cubic  curve  ( BCc  in  short)  defined  on  [0, 1]  with 
its  control  polygon  {Po,Pi,P2,Ps}  and  representation  (hith,2',<pi,<p2',l)-  For 
t € [0, 1],  the  de  Casteljau  algorithm  gives  the  points  q’ ■ defined 

by  the  relation  pjf+ ^ = (1  — t)Pjf  ^ + tPk+k  and  we  obtain  P(t)  = Pq3\  Now, 
we  set  Vfc  = 1, 2, 


7fcJ)  = (Pk-iPk-vPk-iPk-i)  - 


«il) 


j(i) 


Pk+uP&PP), 


and 


P^P-P* 

pWpW 


Ai°>  = 


I pt1)  p(2) 

I rk-irk-i 

I Pk-ilfh 


v>(ix)  = (p1(1)p2(1),p0(1)p1(1))  , 


p(2)  p(3) 

ro  -'  o 

p(J)  p( 2) 

20  1 o 


7i 


(2) 


= (p0(2)p0(3), p0(1)p0(2)) , s‘2)  = (Pj(2)P2(1), P0(3)P<2))  . 


After  some  calculations,  we  obtain  Vfc  = 1, 2, 

= ((1  - t)2  + 2f(l  - t)hk  cos<pk  + ( hkt )2) 2 . 

Then, 

exp(i7fc1))  = ((1  ~t)  + hkzxp{iipk)t) /A{°\ 
exp  (iSk])  = ((!  - f)exp(i<Pfc)  + hkt) /A{°\ 

= 721}  + ^i1}  [27rl- 


Ap  = ((1  - t )2  + 2t(l  - t)/i^  cosy?^  + (fc^t)2)  2 , 


(1) 

(2) 

(3) 

(4) 

(5) 


(6) 
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exp(i7f : >)  = ((1  -t)  + fc^exp^1^)  /A^,  (7) 

exp(*6p^)  = ((1  - t)exp(i^^)  + h^tj  /A^\  (8) 


Remark  1.  For  t € [0,1],  the  de  Casteljau  Algorithm  splits  the  curve  P into 
two  segments:  the  left  one  is  the  restriction  of  P to  the  interval  [0,  t]  defined 
by  the  control  polygon  {Po  , PqX\  Pq2\  Pq3>}  which  admits  the  representa- 
tion (A[0\  A^;7j1^,7p^;it),  whereas  the  right  one  is  the  restriction  to  [t,  1] 
defined  by  {Pq3\  Pj2\  P^\  -P3}  which  admits  {h^  / A^\  h^/A^',  b^^b^', 
A^Aj1^(l  — t))  as  a representation. 


§3.  Curvature  Characterization 

We  determine  the  curvature  radius  R and  the  curvature  p = ^ at  every  point 
of  the  curve  P via  the  parameters  mentioned  above.  Let  us  remember  that 
the  curvature  radius  R is  determined  by  the  approximation  R ~ where 
s denotes  the  arc  length  and  a the  angle  of  the  tangent  vector  with  a fixed 
direction. 

First,  we  give  the  expression  of  R at  t = 0.  Let  us  consider  a value  t near 
0.  We  have  the  following  approximations 


As  ~ (1  + A(x0)  + A^A^Jft, 

Aa  = —7^  — 7j2)  ci  — sin(7j1^  + 7^). 

By  using  (1)  — (7),  we  obtain  when  t tends  to  0 


m = - 


3 

2 


_i 

hi  sint^i 


(9) 


Now  for  t € (0,1],  we  apply  (9)  to  the  right  segment  whose  origin  is  Pq^ . 
Then,  the  curvature  radius  at  t is 


Therefore,  (8)  implies 


m = 


3 A^A^ljl-t) 
2 (h^/Aj1^)  sin«5p^ 


m = - 


3Af\A^)3l 

2 h^  sinipj1^ 


(10) 


Application.  Let  us  consider  the  case  hi  = h%  = h and  p\  = P2  = P with 
p € [0,  f].  One  can  prove  that  Vt  € [0, 1],  A^  = Aj^(=  A),  p^  = p and 
hj1)  = h.  Then,  we  obtain  with  (10),  Vt  € [0, 1],  p{t)  = 2hsin^>/(3ZA4). 

Differentiating  A = ((1  — t )2  + 2t(l  — t)h  cos  p + (ht)2) 2 , we  deduce  that 
hcosp  > 1 and  h < cosiys  are  respectively  the  condition  of  decrease  and 
increase  of  the  curvature  [3]. 


102 


J-C.  Fiorot  and  L.  Schiavon 


§4.  Decreasing  Curvature  Condition  in  the  Case  r*  > 1 

Let  P be  a BCc  curve  with  representation  and  let 

2)  € [0,  |)2.  For  t G [0,1]  , we  set  r = A^/A*0^.  At  t = 1,  this 
parameter  is  r* . Moreover,  we  define 

_ sin(7l(1)  + 7l(2))  _ sin(7j1)  + 7l(2)) 

sin  2^  ’ ^ sin27<2>  ' 

At  t = 1,  we  prove  via  (1)  — (7)  that  these  parameters  are  respectively 


y _ sin(v?i  + <p2) 

sin  2ipi 


sin(ipi  + 1 p2) 
sin  2if2 


Lemma  1.  At  t = 0 and  t = 1,  we  have  respectively  the  equivalences  : 

^ hi  cos  pi  ^1  - T~Y~Sj  > 1 <=>  p'(0)  < 0,  (11) 

h2>^(^-7~Sjcos^p2<=^  p'(l)  <0.  (12) 

Proof:  We  calculate  p'(0)  as  lim;_o  {p{t)  - p(0))/<  via  (9)  and  (10).  We 
obtain 

p'{ 0)  = p(0)(hi(h2s\n(ipi  + <p2)  ~ 2siny?i  - 6sin</Pi(hi  cos</?i  — 1)).  (13) 

Then,  with  the  above  definitions  of  A*  and  p" , we  obtain  (11).  At  t — 1,  we 
consider  the  curve  P(  1 — t)  which  has  the  representation  (l/h2,  l/hi]<p2,<pi; 
h\h2l).  n 


Remark  2.  The  equivalences  (11)  and  (12)  do  not  depend  on  the  interval  of 
definition. 


Lemma  2.  Let  P be  a BCc  curve  defined  on  interval  [i  1 , #2]  with  the  repre- 
sentation {hi,h2',pi,tp2\l)-  Let  us  suppose  that  r*  > 1.  Then, 


p\t  1)  < 0 =»  p'{t2)  < 0. 


Proof:  We  have  the  successive  inequalities 


h2>h\> 


2 

3 


cos  p2  ■ 


The  first  one  uses  the  definition  of  r* . The  second  one  uses  (11).  From  the 
identity 

sin2(</?i  + ip2)  - sin2(<^i  - ip2)  = sin2(/q  sin  2 tp2, 

we  deduce  A * p*  > 1.  After  some  calculations,  we  obtain  the  third  one  and 
consequently  (12).  □ 
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Fig.  1.  Decreasing  curvature  domains  for  r*  = 1, 1.5  and  3. 

Lemma  3.  Let  P be  a BCc  curve  defined  on  [0, 1].  With  the  previous  nota- 
tion and  hypothesis,  we  have 

r*  cos(ipi  + tp2)  — cos 2tpi  > 0 VI  G [0, 1],  r > 1. 

Proof:  For  t € [0, 1],  we  find  that  (A^  A^)2  — (A^)4  is  positive  iff 

(1  - t)2A0  + 2t(l  - t)Ai  + t2A2  > 0 

with 

{Aq  = 2 (r*  cos(y>i  + tp2)  - cos 2 ipx), 

A\  = 2hi(r*  cos ip2  - cosy?i), 

A*  = h\(r*2  - 1), 

is  positive.  Then  r*  cos(y>i  + ip2)  — cos  2^1  > 0.  Conversely,  this  inequality  is 
equivalent  to 

(r*  cos  ip2  — cos  tpi ) cos  ip\  > (r*  sin  <p2  — sin  ipi ) sin  ipi . 

Considering  the  cases  ipi  < <p2  and  ip2  < <pi,  we  see  that  the  coefficients  of 
the  above  Bernstein  polynomial  are  positive.  D 

Proposition  1.  Let  P be  a BCc  curve  defined  on  [0, 1]  with  the  representa- 
tion (hi,  h2;  ipi,  tp2\  l)-  Let  us  suppose  that  r*  > 1.  We  consider  the  domain 

7T 

D i = {((fii,<p2)  € [0,  — ),r*  cos(<^i  + (p2)  - cos  2ipi  > 0}. 

For  (ip\,  tp2)  € D\,  we  have 

p'( 0)  < 0 ^ Vt  € [0, 1],  p'(t)  < 0. 

Proof:  The  proof  is  a consequence  of  Lemmas  1 and  2.  □ 

The  domains  corresponding  to  different  values  of  r*  are  described  via  the 
following  graphs  in  function  of  ip\  (horizontal  axis)  and  ip2  (vertical  axis).  In 
Figure  1,  the  white  part  represents  the  decreasing  curvature  domain  whereas 
the  dark  one  denotes  a domain  where  the  decrease  at  t = 0 is  not  possible 
(r* A*  > 3).  One  can  notice  that  the  latter  is  empty  when  r*  > 6.  In  the  grey 
part,  we  cannot  say  anything  about  the  monotonicity. 
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§5.  Decreasing  Curvature  Condition  in  the  Case  r*  < 1 

With  r*  < 1,  the  condition  p'(0)  < 0 is  too  strict,  so  we  consider  the  following 
sufficient  condition  on  decreasing  curvature  at  the  origin: 

p'(0)<  4(1  — ^)  p(0), 

which  is  equivalent  to 

^ /ucos<pi  ^1  - > r*_1.  (14) 


If  we  set 


(14)  becomes  r*a* 


\ hlCOStfi!  fl  - 


> 1.  Furthermore,  we  set  Vf  € [0, 1], 


3 A(o)  (1)  (,  r\\ 

a = - A\  ; cos  7)  M 1-  y ) . 


A calculation  gives  a = (1  — t)  + a*t. 

Remark  3 . The  inequality  (14)  does  not  depend  on  the  interval  of  definition. 

Lemma  4.  Let  P be  a BCc  curve  defined  on  interval  \ti,t2}  with  the  repre- 
sentation (hi,h2',ipi,tp2',l)  and  r*  < 1.  Let  us  consider  domain 

7 r 

D2  = {(<P1,<P2)  G [0, -[,r*cos<p2  - cost/?!  > 0). 


Then  V (<pi,<p2)  G D2, 

p'{t  1)  < 4 ^1  - p{h)  =>  p'{t2)  < 0 
Proof:  For  (<pi,<p2)  G D2,  we  have 


which,  with  (14),  implies 


t *L  .2  ^ r*\*  1 ^ 3 fl  r> 

h2  = r hi  > - 1 — > ~ 1 1 

3 \ 3 / cos<pi  2 


> A ( 1 3-  ) cos ^2. 


□ 
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Lemma  5.  Let  P be  a BCc  curve  defined  on  [0,1].  Under  the  previous 
hypothesis,  we  have 

V(<pi,tp2)  € D2  &\/t  e [0,1],  rcos7^  -cos  7^  > 0. 


Proof:  For  t € [0, 1],  we  find  that  r cos  — cos  7^  > 0 if  and  only  if 

(1  — t)(r*  cos(y>j  + ip2)  — cos2<pi)  + h2(r*  cos  ip2  — cos<pi)t  > 0. 

Then  r*  cos  <p2  — cos  > 0.  Conversely,  r*  cos  ip2  — cosipi  > 0 and  r*  < 1 
imply  <p2  < ip\.  Then  r*  sin  <p2  — sin  <p2  < 0.  Consequently, 

r*cos(ipi-Up2)— cos  2ipi=(r*  cos(<pri(p2)— cos2ipi)cos<pi— (r*smip2— sin(/3i)sin^i 

is  positive.  □ 

Lemma  6.  Let  P be  a BCc  curve  defined  on  [0,1].  Under  the  previous 
hypothesis,  we  have  V(<p\,ip 2)  € D2, 

hi(r*2a*  — 1)  + 2 (r*a*2  cosip2  — cosipi)  > 0 =>  Vt  € [0, 1],  ra  > 1. 

Proof:  The  last  inequality  is  equivalent  to  deciding  whether  a fifth- degree 
polynomial  is  positive.  We  verify  that  all  its  coefficients  are  positive  but  for 
one.  The  positivity  of  this  coefficient  is  equivalent  to  the  first  inequality  in 
the  lemma.  □ 

Proposition  2.  Let  P be  a BCc  curve  defined  on  [0, 1]  with  the  representa- 
tion (hi,h2;<pi,<p2;l).  Let  us  suppose  that  r*  < 1.  We  consider  the  domain 
D2  as  mentioned  above  and 


D3  = € [0,  | ),hi(r*2a * - 1)  + 2(r* a*2  cos  ip2  - cos  ipi)  > 0}. 

Then,  V(ipi,ip2)  € D2  n D3, 

p'( 0)  < 4 (l  - p(0)  =>  Vf  € [0, 1],  p'(t)  < 0. 

Proof:  The  proof  is  a consequence  of  Lemmas  5,6  and  4.  □ 

Here,  we  describe  the  different  admissibility  domains  D2f)D3  (represented 
in  white)  for  several  values  of  r*  in  considering  that  a*  = 1 /r* . As  illustrated 
in  Figure  2,  there  is  a domain  continuity  when  r*  is  near  1 with  the  case 
r*  = 1 (Figure  1).  When  r*  decreases,  the  domain  gets  smaller  and  smaller 
and  then  it  is  finally  empty  for  r*  ~ 0,33.  If  we  take  a*  1/r* , the  domain 
size  increases  and  tends  to  D2. 
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Fig.  2.  Decreasing  curvature  domains  for  r*  = 0.99,0.95,0.9  and  a*  = 1 /r*. 

§6.  Examples 

Our  results  are  illustrated  by  examples  of  BCc  curves  whose  curvature  de- 
creases. For  each  curve,  one  quarter  of  the  curvature  radius  is  represented  at 
t values  0,  0.1,  0.2,  ...,  1. 

Example  1.  r*  = 1 and  p'(0)  = 0.  The  first  curve  with  small  angles  (f , 
and  a length  rate  hi  ~ 0.99  (Figure  3 (a))  is  characterized  by  a small  radius 
increase.  The  second  one  has  angles  ( |)  and  hi  ~ 1.36  (Figure  3 (b)). 

Example  2.  (<pi,<P2)  = (§,f)  and  p'(ty  = 0.  The  curve  with  r*  = 1.25 
(Figure  4 (b))  is  tighter  than  the  curve  with  r*  = 1 (Figure  4 (a)). 

Example  3.  r*  = 1.5  and  (<pi,<p2)  = We  compare  the  curvature 

radius  increase  at  the  origin  for  a curve  (Figure  5(a))  with  p'(0)  = 0 (hi  ~ 
1.71)  and  another  one  (Figure  5 (b))  with  p'( 0)  = —0.12  p{ 0)  (hi  ~ 2.05). 

Example  4.  We  consider  two  curves  with  r*  < 1 and  a*  = 1 jr*\  r*  = 
0.95,  (^1,^2)  = (f,f)  and  r*  = 0.9,  (<^1,^2)  = The  parameters 

obtained  by  calculations  (take  the  equality  in  (14))  are  respectively  hi  ~ 1.44, 
p'(0)  = —0.03  p(0)  and  hi  ~ 1.28,  p'(0)  = —0.06  p(0)  (Figure  6 (a)-(b)). 


Fig.  3.  (a)  r*  = 1,  (<pi,<p2)  = (f,  fe)  > (b)  r*  = b (VUV 2)  = (f , f )• 
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